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Abstract 

We consider stochastic differential equation involving pathwise integral with re- 
spect to fractional Brownian motion. The estimates for the Hurst parameter 
are constructed according to first- and second-order quadratic variations of ob- 
served values of the solution. The rate of convergence of these estimates to the 
true value of a parameter is established. 
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1. Introduction 

Consider stochastic differential equation 

Xt^^+ f f{X,)ds+ f g{Xs)dBf, te[o,T], r>o, (1) 

where / and g are measurable functions, is a fractional Brownian motion 
(fBm) with Hurst index l/2<ff<l,^isa random variable. It is well-known 
that almost all sample paths of have bounded p- variations for p > \/H. 
Therefore it is natural to define the integral with respect to fractional Brownian 
motion as pathwise Riemann-Stieltjes integral (see, e.g., [l6| for the original 
definition and for the advanced results). 

A solution of stochastic differential equation (jlj on a given filtered probabil- 
ity space (17, J',P,F = {Ft},t& [0,T]), with respect to the fixed fBm (B^,F), 
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1/2 < H < 1 and with J^j-measurable initial condition ^ is an adapted to the 
fihration F continuous process X = {Xt : ^ t ^ T} such that Xq = ^ a.s., 

\f{Xs)\ (^^+\ '^■Sf I < oo^ = 1 for every Os^ts^T, 

and its almost all sample paths satisfy ([1]). 

For < a < 1, Ci+"(R) denotes the set of all C^-functions 5: M M such 
that 

sup + sup ^^^^-^^ < K). 

X x^y \X y\ 

Let / be a Lipschitz function and let g G C^+"(M), - 1 < a ^ 1. Then 
there exists a unique solution of equation ([T]) with almost all sample paths in 
the class of all continuous functions defined on [0, T] with bounded p-variation 
for any p > jj^ (see [13], [l3] and [10]). Different (but similar in many fea- 
tures) approach to the integration with respect to fractional Brownian motion 
based on the integration in Beso v sp aces and corresponding stochastic differ- 
ential equations were studied in 15|, see also [a] and [14] where the different 



approaches to stochastic integration and to stochastic differential equations in- 
volving fractional Brownian motion are summarized. 

The main goal of the present paper is to establish the rate of convergence 
of two estimates of Hurst parameter to the true value of a parameter. The 
estimates are based on the two types of the quadratic variations of the observed 
solution to stochastic differential equation involving the integral with respect 
to fractional Brownian motion and considered on the fixed interval [0,T]. The 
paper is organized as follows: Section 2 contains some preliminary information. 
More precisely, subsection 2.1 describes the properties of p- variations and of the 
integrals with respect to the functions of bounded p- variations while subsec- 
tion 2.2 contains the results on the asymptotic behavior of the normalized first- 
and second-order quadratic variations of fractional Brownian motion. Section 3 
describes the rate of convergence of the first- and second-order quadratic varia- 
tions of the solution to stochastic differential equation involving fBm. Section 4 
contains the main result concerning the rate of convergence of the constructed 
estimates of Hurst index to its true value when the diameter of partitions of the 
interval [0, T] tends to zero. Section 5 contains simulation results. 



2. Preliminaries 

2.1. The functions of bounded p-variation 

First, we mention some information concerning p-variation and the functions 
of bounded p-variation. It is containing, e.g., in Q and [l^. Let interval 
[a, b] C M. Consider the following class of functions: 

Wp{[a,b]) := {/ : [a, 6] ^ M: Vp{f; [a,b]) < 00}, 
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where 



Vpif; [a,b]) = sup^\f{xk) - /(a;fe-i)|^. 



fe=i 



Here tt = {xf. i = 0, . . . ,n} stands for any finite partition of [a,b] such that 
a = xo < Xi < ■ ■ ■ < Xn = b. Denote n([a, b]) the class of such partitions. We 
say that function / has bounded p-variation on [a; b] if Vp(f: [a, b]) < oo. 



LetVpif) ■.= Vp{f;[a,b]) 



,1/p 



(/; [a, b]). Then for any fixed / we have that 



Vp{f) is a non-increasing function of p. It means that for any < q < p the 

relation Vp{f) ^ Vq{f) holds. 

Let a < c <b and let / G Wp([a, 6]) for some p € (0, oo). Then 

Vp{f; [a,c]) +Vp{f; [c,b]) < Vp{f; [a,b]), 
Vp{f; [a,b]) ^ Vp{f; [a,c]) + Vp{f; [c,b]). 

Let / G Wq{[a,b]) and h G Wp{[a,b]), where p''^ + q'^ > 1. Then the 
well-known Love- Young inequality states that 



f fdh-f{y)[h{b)-h{a)] 

J a 



whence 



(2) 



(3) 



Here Vq^oo{f][a,b\) = [«- ^]) + sup^^^^b Cp,g = C(p ^ + q ^) and 

C(s) = X^n^i is the Riemann zeta function. Further, for any y G [a, b] 

Vp(^J^ [f{x)-fiy)]dh{xy,[a,b]^ ^Cp,g[Vg{f;[a,b]) 

+ sup \f{x)-f{y)\]Vp{h;[a,b])^2Cp,gVg{f;[a,b])Vp{h;[a,b]). (4) 



Denote 



1^1 



oo ^ sup|A(a;)|, = sup 



la; -2/1" 



Let F be a Lipschitz function and let G G Ci+"(ffi) with < a < 1 and 
l^p<l + a. Then for any h G Wp([a, 6]) 

Vp,^{F{hy, [a,b]) ^LVp{h; [a,b]) + sup \Fih{x)) - F{h{a))\ + \F{h{a))\ 

^2LVp{h:[a,b]) + \Fih{a))l (5) 
Vp/a,oo{G{h); [a,b]) ^Vp,oo{G{h); [a,b]) ^ 2\G'\^Vp{h; [a,b]) + \G{h{a))\, (6) 
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and 



^p/a,oo(G'W;[a,6]) 

< [a,b]) + sup \G'{h{x)) - G'{h{a))\ + \G'{h{a))\ 



<2|G'UV;,"(k [a,b]) + \G'{h{a))\. 



(7) 



Let / e Wp([a, 6]) and pi > p > 0. Then 



(/; [a, b]) < Osc(/; [a, 5]), 



(8) 



where Osc(/; [a,b]) = sup{|/(a;) - f{y)\: x,y £ [a,&]}. 

Take functions /i,/2 G Wp{[a,b]), < p < 00. Then /1/2 € )%([a,6]) and 



Vp{fif2; [a,b]) < yp,oo(/i/2; [a, 6]) < CpVp,^{fv, [a,b])Vp,oo{f2; [a,b]). (9) 



Let /i G Wq([a, 6]) and /2 G Wp([a, 6]). Then it follows from Young's version of 
Holder's inequality that for any partition tt G n([a, 6]) and for any p~^+q~^ ^ 1 



Y,Vg{h;[xi-uXi])Vp{f2;[xi-uXi]) < Vg{fv,[a,b])Vp{f2;[a,b]). (10) 



Second, we state some facts from the theory of Riemann-Stieltjes integration. 
Let / G Wq{[a,b]) and h G Wp{[a,b]) with < p < 00, q > 0, 1/p+l/q > 1. Let 
symbol (R) stands for the Riemann integration, and (RS) stands for Riemann- 
Sticltjcs integration. Then integral (RS) J^^ f dh exists under the additional 
assumption that / and h have no common discontinuities. 

Proposition 1. Let f : [a,b] M. be such function that for som,e 1 ^ p < 2 
f G CWp([a, 6]). Also, fei F : M — > M 6e a differentiable function with locally 
Lipschitz derivative F' . Then composition F'{f) is Riemann-Stieltjes integrable 
with respect to f and 



Furthermore, the following substitution rule holds. 

Proposition 2. Let /i,/2 and fs be functions from CWp{[a,b]), 1 ^ p < 2. 
Then 



{RS) / f,{x)di{RS) / f2{y)dfs{y)]={RS) / fi{x)f2{x)dfi{x). 



F{f{b)) - F{f{a)) = (RS) / F'ifix)) df{x). 







Finally, assume that 




and 
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where /i is continuous function, /2,/3 G CWp{[a,b]) for some 1 ^ p < 2, and 
Q is a differentiable function with locahy Lipschitz derivative q. It follows from 
Propositions [T] and [2] that 



Q{F,{x) + F^ix)) - Q{0) ^ / qiF,iy) + F2iy))diF,iy) + F2iy)) 

J a 

q[Fi{y) + F-2{y))h{v)dy+ I q[Fi{y) + F2{y))f2[v) dh{y). (11) 



2.2. Asymptotic property of the first- and second-order quadratic variations of 
fractional Brownian motion 
Consider the fractional Brownian motion (fBm) = {B^ ,t £ [0,T]} with 
Hurst index H S (^j !)• Its sample paths are almost all locally Holder up to 

order H. Moreover, for any < 7 < we have that Lj,'^ :— sup s^t H-sh 

is finite a.s. and even more, E(L^''*^) < 00 for any k ^ 1. The following 
estimate for the p- variation of fBm is evident: 

^i^'^/^(t-s)i/p, (12) 

where s < t ^ T, p > 1/H. 

Let 7r„ = {0 = to < ^1 < ' ' ' < *n = T}, T > 0, be a sequence of uniform 
partitions of interval [0, T] with t^ = ^ for all n G N and ah fc e {0, . . . , n}, 
and let X be some real- valued stochastic process defined on the interval [0,r]. 

Definition 3. The normalized first- and second-order quadratic variations of X 
taking along the partitions (7r„)„gN and corresponding to the value 1/2 < H < I 
are defined as 



a(1) 

k=l 

and 



ri-l 
k=l 

For simplicity, we shall omit index n for points t^ of partitions 7r„. 

It is known (see, e.g., Gladyshev JJ) that Vn^\B^ ,2) — >■ T a.s. as 71 00. 
Also, it was proved in Benasi et al. Q and Istas et al. [§] that Vn'^\B^ ,2) — >• 
(4 - 2'^")T a.s. as n ^ 00. 

Denote 

r(t) 

,(1) ^H^2 



y«(i?^,2), =n^^-i^(AWB«)= 

k=l 

yp)(i?^,2), =n^^-l'|]\AgB^)^ t e [0,T], 



fe=i 



5 



where r{t) — max{fc: tk < t} — [^]. It is evident that 

BV^'Hb" ,2)t ^ Pit) and EVJi'\b'' ,2)t ^ (A - 2''')p{t), (13) 

where p{t) = m.a.x{tk : ^ t}. 

The fohowing classical result will be used in the proof of Theorem [Sj 

Lemma 4. (Levy-Octaviani inequality) Let Xi, . . . , X„ be independent random 
variables. Then for any fixed t, s ^ 



P I max 



>t + s\ 



P( 




>t) 


1 - maxi^i^„ P( 


>s) 



□ 



Theorem 5. The following asymptotic property holds for the first- and second- 
order quadratic variations of fractional Brownian motion: 

sup 2)t -EV;(^)(B^,2)J = 0(71-1/2 In^/^^) i = l,2. (14) 

Remark 6. It follows from ^ and ^ that sup^^j. \vJi'\B" ,2)t\ = 0(1) 
a.s., i — 1,2.. 

Proof. We can consider yi*-* (^B^ , 2) as the square of the Euclidean norm of 
the n-dimensional Gaussian vector X„ with the components 



1 s^k s^n- (i-l). 



Obviously, one can get a new n-dimensional Gaussian vector Xn with inde- 
pendent components applying the linear transformation to Xn. It means that 
there exist nonnegative real numbers (A^*^„, . . . , A^'^^^_j^j ^) and such n — (i — 1)- 
dimensional Gaussian vector Yn with independent Gaussian A/'(0, l)-components 
that 

n-(i-l) 

yW(s«,2)^ Y: A,,„(ri^))^ 

The numbers (A^'|j, . . . , n) eigenvalues of the symmetric n— [i — 

1) X n — (i — l)-matrix 

\ ' / l^j,/c^n— (2 — 1) 

Now we can apply the Hanson and Wright's inequality (see Hanson et al. 
or Begyn [1]), and it yields that for e > 



Ji-(i-l) 

j=k 



> e ) ^2 exp ^— min 
^2 exp — min 



, *(j) ' Y^n-(i-l) /, (i) 



Cie 



(15) 
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where Ci, C2 are nonnegcitive constants, ^ — ^^^k^j^n—{i~i) — 
The evident equahty holds: 

n-(i-l) 

Furthermore, it follows from that the sequence EVi*-* (S^, 2), n > 1 is 

bounded. So, the sums X]j'=i*~^' '^j'i ^'''^ bounded as well. It is easy to check 
that 

n— (i — 1) n — (i— 1) 

Therefore for any < £ 1 the inequality ([T5|) can be rewritten as 

n-(i-l) 



> e ) ^ 2 exp 



a; 



(16) 



where is a positive constant. 

Now we use Levy-Octaviani inequality (see LemmaH]) and evident inequality 



l-x 



^ 2a; for < a; 1/2 



to obtain the bound 
P 



max 

l^/£^n-(i-l) 



> 2e 



^°p(-g^) 



1 — 2 exp 



^ 4 exp 



assuming that 



exp 



An 



< 1/4 and < e < 1. 



So, for the values of parameters mentioned above. 



P\n'^"~^ max 

\ l<fc<n-(i-l) 

Furthermore, 



E(A<:is^)^-EE(A<:ii?^)^ 

i=i j=i 



> 2e ^4 exp 



n-(i-l) 



a: 



;«<AV«-i max E Mini' 

" l<k<n-(i-l) ^ ' 
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where = l\f.^B" . From Gladyshev fzll and Begyn [H we get 

a;(') < C7r\ (17) 

Now we set 



and conclude that 



max 

l^fe^n-(i-l) 



2 2^ —1 1 

e„ = — — n mn 



4 



> 2e, 



^ 4 exp —2 Inn 



It means that 

oo / 



2ff-l 

n max 



l^fc^ri-(i-l) 



j=i i=i 



> 2e„ < cx). 



Finally, we get the statement of the present theorem from the Borel-Cantelli 
lemma and the evident equality 



sup|yW(i?^,2)t-Eyr(i3",2) 



2H~1 

n max 

l^fc^n-(i-l) 



3. The rate of convergence of the first- and second-order quadratic 
variations of the solution of stochastic differential equation 

First, we formulate the following result from fTJ] about convergence of first- 
and second-order quadratic variation. 

Theorem 7. Consider stochastic differential equation ([T]), where function f is 
Lipschitz and g G C^+" for some < a < 1. Let X be its solution. Then 

fT 

(X, 2) cW / g^{X{t))dt a.s. as n oo, (18) 

where 

cW^/l for 1 = 1, 

|(4_22ff) fori = 2. 

Second, we prove the following auxiliary result. 



8 



Lemma 8. Let X he a solution of stochastic differential equation ([T]). Define 
a step-wise process X'^ that is a discretization of process X: 



X^ ^ 



X{tk) for te[tk,tk+i), 0,1,..., n- 2, 

X{tn-l) for te[tn~l,tn]. 



Then for any p > ^ we have that 



svip\X[ -Xt\ = 0{n-'/P). 

t<T 



Proof. Consider t G [tk,tk+i). We get immediately from the Love- Young 
inequality 1^ that 



1X7 - Xt 



f{Xs)ds 



P"(i) 



9(Xs)dB^ 



H 



^Tn-^ sup \fiXt)\+Cp,pVp,MX);[tk,tk+i])Vp{B";[tk,tk+i]). 



Further, 

sup \f{Xt)\^snp\f{Xt)~f{O\ + \fm^LVp{X;[0,t]) + \f{Ol (19) 

where L is a Lipschitz constant for /, and 

Vp,oo{9{X);[tk,tk+i]) ^\g'\ooVp{X;[tk,tk+i]) + sup\g{Xt)^g{0\ + \9{0\ 



i^2\g'\^Vp{X;[0,T]) + \g{O]. 



(20) 



We get the statement of the lemma from (IT^ , (PH)) and inequality (IT^ . ■ 

Now we prove the main result of this section which specifies the rate of 
convergence in Theorem [T] 

Theorem 9. Let the conditions of Theorem^hold and, in addition, a > jj—l- 
Then 



y«(X,2)-c« / g^{X{t))dt = 0{n-^/Hn'^^n). 



(21) 



Proof. Decompose the left-hand side of (pT|) into three parts: 



/« :=Vj^iX, 2) - c« / g^Xit)) dt = L^^'^^ + L^"^ + 
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where 

n—{i—l) n—{i—l) 

fc=l k=l 

n-(i-l) 

fc=i 

ri-(i-l) 

- y: 5^(^.-l+(.-l))E(A«s^)^ 
fc=i 

ri-(i-l) rp 

k=l -^0 

and Xfc = X{tk)- We start with the most simple term In'^'' and get immediately, 
similarly to bounds contained in ^ , that for any P > jj 

fc = l-(i-l) ''tk-l+(^-l) 

sc:2c«r|5'U sup - Xt\ ■ [WlooVpiX; [0, T]) + |g(0|] . 
In order to estimate /n^'*'' , denote 

rW-(i-l) 

fe=i 



Then 



n-(i-l) 

k=l ^0 



and 



71 



n-(i-l) 

^"-^ E 5^(^.-i+(.-i))[(Ai:>-)^-E(Ai:>-)^] 

T 
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Note that 1/p + 1/2 > 1 ^ < p < 2. Therefore, we obtain from the 
Love- Young incquahty and from ([H])-® that 



Jo 



^Cp,2Vp,^{g^{X); [0,T])l/2(^(^' - E5«; [0,T]) 
s=:Cp^2{Osc(5« - E5«; [0,r])}'/Vp^oo(<?'W; [0,T]) 
X -E5W;[0,r]) 

^2CpJsup - E5«|)'^Vp2^(g(X); [0,r]) 



,2//- 



1/2 



It follows from Theorem [5l Remark IH and ([6|) that the rate of convergence 

It still remains to estimate In ' ■ Consider only i = 2, the proof for i = 1 is 
similar. 
Denote 



4- 



- 



[/(X,) - J{Xk)] ds 



tk 
tk 



tk 



[f{X,) - f{Xk)] ds, 



tk-l 



tk- 
ik + 1 

tk 



tk 



g{Xk)-g{X,)^ g'{Xk)f{Xk)du- 



j\'{Xk)g{Xu)dB^^dBf, 



g(Xs)~g{Xk)- g\Xk)f{Xk)du~ I g' {Xk)g{Xk) dB'^ ]dB, 
Jtk 

tk4-^ ftk 



tk 



s ' 



4 ^g'{Xk)f{Xk) ( / (s - tk) dB^ + (tk- s) dBf) 



- 



\g'{Xk)g{Xk){{I^^^lB"f + { 



> B^ 



(2) 



Equalities 

Ctk 



tk-1 
tk-\-l 



tk 



tk 



g'{Xu)g{Xk)dB^ ]dB^ =i ,g'(Xfc)ff(Xfe)(A[?'i3^)^ 
g'iXk)giXk)dB^ )dB^ =i g'(X,).g(Xfe)(A«,_„S^)l 



4^giXk)A^^>B 



(see Proposition [T]) imply 



11 



Taking into account Lipschitz property of / and Lemma[51 we can conclude that 
\f{Xt) - f{Xk)\ s=: Lsup ^Xt\^ 0(n-i/P). 



Therefore 



n-1 „t 



fe=l 



[fiXs)~fiXk)]^ds 



fc=i 

rtk 



fc = l "'*'=-! 



Consider J|. It follows from equality ((TT|) that for any fixed t G [t/c-i, ifc] 
giXk) - g{Xt) = I ' g'{X,)f{X,) ds + T g'{X,)g{X,) dBf . (22) 



Substituting equality ([^^ into J| we get 



tk 



tfc-1 •'S 



g'{X^)f{Xu) - g'{Xk)f{Xk)]dudB^ 



(23) 



Transforming identically the first term in the right-hand side of and 
applying to it Love- Young inequality ([2]) , we conclude that for any p > jj 



g'{Xu)f{Xu) - g'{Xk)f{Xk)]dudBf 
^Cp^iVi(^J ' [g'{Xu)fiX^) ^ g'{Xk)f{Xk)]du; [tk-i,tk]^Vp{B" :[tk-i,h]) 

rtk 

^Cp,iVp{B'';[tk-i,tk]) \g'{X^)f{X^)-g'{X,)f{Xk)\du. (24) 



Henceforth we consider the following interval of the values of p: < p < 1 -f- a. 
Then it follows from inequality (j4|) that the second term in the right-hand side 
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of admits the bound: 



[g'(Xu)9{X^) - g'{Xk)g{Xk)]dB^ dB^" 



tk-l ■J s 

^Cp^p/cVp/J I [g'{X^)g{X^) - g'{Xk)g{Xk)]dB^;[tk-i,tk[ 



X Vp{B";[tk-i,tk]) 
< 2Clp^^Vp/4g'{X^)g{X^y, [tk^,,tk])Vp' {b" ; [tk-i,tk]). (25) 



We conclude from (I23])"(125]) that 

|J|| ^Tn-^Cp,,Vp^4g'{X)f{X); [tk,tk+i])Vp{B" ; [tk,tk+i]) 
+ 2Clp/^Vp/4g'{X)g{X);[tk,tk+i])Vp^{B";[tk,tk+i]). 

Applying inequalities ^ and (ITUl) we immediately obtain that 



^2T^Cl,n-\ ma^ [Vpf4g'{X)fiX);[tk,tk+i])Vp{B'';[tk,tk+i])] 



k=l 

X Vp/4g'iX)fiX); [0,T])Vp{B"; [0,T]) 

+ 4CpV ^a'^ , [Vp/4g'iX)giX);[t,,tk+i])Vp\B";[tk,h+i])] 

X Vp^4g'{X)giX); [0,T])Vp{B" ; [0,T]) 
^2T'Cl,n-\max [Vp{B";[tk,tk+i])]Vp%^^{g\Xy,[0,T]) 

X v;2^(/W;[o,T])yp(s«;[o,T]) 

X [0, T])Vp%^^ {g{X); [0, r])yp(i3^; [0, T]). 

It follows from the inequalities ©-([I]) that the values of the variations 
yp,oo(/(X);[0,r]), Vp^^^^{g{X);[0,T]), and Vp/^^^{g'{Xy,[0,T]) 
are finite. Therefore we get from that 

n 

fe=i 

The similar reasonings lead to the similar bound for J|, and we conclude 
that 

k=0 

Consider J^. It consists of two terms that can be estimated in a similar way. 
Applying inequalities ([2]) and ([T2l) . we obtain the following bound for the first 
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term: 



fe=i ^ 

ra-l 



k=l 



< n-^T^C^i max [5'(X,)/(Xfe)yp(B^; = ©(n-^-^/f). 

As a consequence, 

n-l 

fe=l 

Furthermore, note that under our assumptions sup^gjQ j^] < cxd and 

supgg[o,T] l5'(-'^s)l < 00 a.s. Therefore we have for the first term in J| that 

n—\ n— 1 

k=l ^ ^" /c=0 

The second term is bounded in a similar way, and we conclude that 

n-l 

E(J|)' = oK-^/^)- 
fe=i 



Thus 

5 N 2 



(Y^A =0{n-'-^/PVn-^/PVn'-^/P)=0{n'-^/P). 
^ 1=1 ' 

At last, 

n-l 



fe=l 



=0(n-4/f+2^) (26) 

for any ;i= < p < 1 + a. Set 1/p = - £ for £ < (ff/2 - 1/16) A (if - -^). 
Then 



= 0{n-^l^hi^l^n). (27) 
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4. The rate of convergence of estimators of Hurst index 

Consider the following statistics: 

and construct the following estimate of Hurst index H: 



1 1 

where 

Cn = |2-i(l-2n-i/4(lnn)i/4+/3) <; <; l + 2n-^/\\nn)'/*+py /? > 0. 

Further, introduce the following notation: g^'^\T) — c*^*^ g^{Xs)ds. 

Theorem 10. Let conditions of Theorem^hold with a > jj — 1. Also, let X 
be a solution of (Qp and assume that random variable (T) is separated from 
zero: there exists a constant cq > such that g^'^\T) ^ cq a.s. Then Hn ^ is 
a strongly consistent estimator of the Hurst index H and the following rate of 
convergence holds: 



-H\ = 0{n-^/\lnny/*+^) a.s. 



for any /3 > 0. 



Proof. Consider a sequence 1 > (5„ 4- as n — s- cxd. It will be specified later on. 
Introduce the events 

Cn - ^1(1 - Sn) ^ R^:^ ^ 1+(5„|. 

Also, introduce the notations 

A^^ ^V2^'^\X,2) and B«=y«(X,2) 
and note that 2'^"-^!^^ = Then 

and estimate Hn^ has a form 



'2 2 In 2 " ' ^" 
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It is easy to see that C„ := f2\C„ has a form 



< 2 



2JJ-2 



(1 - 5n) 



> 2 



(1 + (5„; 



c 



(28) 



Then 



^W-ffl ^ ■ (2nr-^y,(:)(X,2) 

2 In 2 



The latter representation implies that 



\H^^ -H\ <H1 



+ 



{|fT-i|>^"} 21n2 



In 



/ 1 



In- 



V21n2 {2--=(i-5„)^^<i-5„} 
1^21n2 ""Bi'^J {i+5.^4t<2'"-Mi+^")} 



E^»- (29) 



1=1 



In what follows we need an elementary inequalities: — ln(l— x) < 2 ln(l+x) < 2x 
provided that 0^x^1/2. 

Consider L^. We divide it in two parts. As to the first part, it is obvious 
that 



In 



A 



(0 



and 



In 



1 



1 



A. 



(i) 



B. 



An' A 

1 — (5« ^ — 7^ < 1 implies that < 1 



B 



{'') 



B 



(i) 
n 



Applying inequality — ln(l — a;) ^ 2x, < a; ^ 1/2, we deduce that for 5n ^ 1/2 



-In- 



A. 



As to the second part, 



^ 2 1 



Aii) 



^ 2<5„1 



{i-.„.^<i}- 



In 



A(^) 



In 



On 
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Consider Lf^. From here we easy deduce that 



1 



In 



A. 



21n2 5 W / {22ff-2(i-5„)s:^<i-5„} 

21n2'- ^ |22-H-2(i-5„)s;^<i-5«| 



^ 1 - 



ln(l - 5,,) 
2 In 2 



"{22H-2(i_5„)^^<l-5„} 



hl2y |22H-2(i-A-„)5;^<i_5„} 

(l + 2<5„)l,^w 

The term L^^ is estimated similarly as the second part of L^. Thus we get 



In 



A. 



A. 



11, 



<(1 + 25„)1,^(„ 



Summarizing, we conclude that 



S=:(1 + 2J„)1,|^(,) I .+2Sn. 

Now, let /3 > 0. Note that 



n 



- 1 



(i) 



> 5n,B^^ ^ (lnn)-^|U{B,« < (Inn)"''} 

= {|A« - B«| > <5„i3«,i3« ^ (Inn)-^} U < (Inn)-'^} 

C - bW| > ^„(lnn)-^} U < (Inn)-^}. 

Therefore 

\Bn^ - H\ ^ + 2^")l{|A<:'-S<;>|><5r.(lnri)-(3}u{B<.''<(lnn)-/3} 

It follows from ^T^i that 

|AW-B«| = 0(n-V4inV4,) 



+ 2(5„ 
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and 



g^{X,)ds 



= 0(n-i/4lni/4n). 



Obviously, for any n > exp{(^) } we have that g'^^\T) ^ cq > 2(lnn) ^ a.s. 
and 

{B« < (Inn)-''} {SW < (lnn)-'3,g«(T) ^ 2(lnn)-'^}. 
Now, let 5n < (Inn)-'^. Then it is not hard to deduce that 

(B^*) < (lnn)-'',gW(r) ^ 2(lnn)-^} 

= < (lnn)-^g«(T) ^ 2(lnn)-^i?« < 5«(T) - 5„} 

C{|i3«-5«(r)|>5„}. 

Therefore, 

< (Inn)-^} C {|B« -5^')(T)| ><5„} 

ifn>exp{(|^)^}. 

Finally, specify (5„. More precisely, set (5„ = n~^/'*(lnn)-'^/^+^'', ^ > 0. Note 
that Sji < (lnn)~^ for sufficiently large n and, moreover. 



<5„(lnn)-^ ~ n-i/4(lnn)i/4+/3 







a.s. as n — > oo. 



The latter relation together with Theorem [9] imply that for any a; G fi' with 
P(ri') = 1 there exists hq = uq^uj) such that for any n > uq 

and we obtain the proof. 



5. Simulation results 

Consider fractional Ornstein-Uhlenbeck process that is the solution of the 
linear stochastic differential equation 

dXt = -Xtdt + dSf , Xq = 0. 

with the step 0.05 and for increasing (in the logarithmic scale) number n of 
points from n = lO'^ to n = 10^. Table 1 presents the values of the difference 
- H\ for the values of H from 0.55 to 0.95. We can conclude that the 
difference \Hn'' — H\ decreases rapidly in n and for fixed value of n increases in 
H . Table 2 demonstrates that the rate of convergence agrees with Theorem [TOl 
at least, for /? = 0.05. Moreover, we can see from Table 3 that in the case of 
the linear equation the rate of convergence for H E (0.5, 0.7) can be estimated 
by n-i/2(lnn)i/2. 
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Table 1: 0P ~ H\ 



H 


100 


250 


1000 


2500 


lO* 


2.5 ■ 10* 


lO-""' 


2.5 ■ lO"' 


10** 


0.55 


0,08401 


0,05488 


0,02124 


0,01467 


0,00777 


0,00549 


0,00195 


0,00160 


0,00079 


0.6 


0,07216 


0,04145 


0,02213 


0,01286 


0,00683 


0,00466 


0,00214 


0,00137 


0,00069 


0.65 


0,07761 


0,04811 


0,01972 


0,01296 


0,00626 


0,00414 


0,00210 


0,00144 


0,00066 


0.7 


0,05364 


0,03403 


0,02023 


0,01219 


0,00608 


0,00341 


0,00183 


0,00125 


0,00065 


0.75 


0,06485 


0,03798 


0,02187 


0,01147 


0,00707 


0,00424 


0,00211 


0,00140 


0,00083 


0.8 


0,05938 


0,03884 


0,02040 


0,01307 


0,00791 


0,00528 


0,00303 


0,00227 


0,00120 


0.85 


0,04666 


0,03577 


0,02105 


0,01684 


0,01011 


0,00753 


0,00511 


0,00384 


0,00249 


0.9 


0,06311 


0,04642 


0,03037 


0,02338 


0,01667 


0,01352 


0,00984 


0,00801 


0,00599 


0.95 


0,06219 


0,04763 


0,03488 


0,02907 


0,02295 


0,02018 


0,01640 


0,01448 


0,01213 



Table 2: 1^4^^ - H\ ■ nO-2S(lnn)-°-3 



n points 



H 


100 


250 


1000 


2500 


10* 


2.5 • 10* 


10^ 


2.5 • 10^ 


10« 


0.55 


0.16802 


0.13070 


0.06689 


0.05596 


0.03991 


0.03449 


0.01663 


0.01681 


0.01136 


0.6 


0.14433 


0.09873 


0.06969 


0.04906 


0.03506 


0.02923 


0.01828 


0.01435 


0.00994 


0.65 


0.15522 


0.11458 


0.06212 


0.04942 


0.03215 


0.02602 


0.01793 


0.01510 


0.00954 


0.7 


0.10727 


0.08104 


0.06370 


0.04652 


0.03125 


0.02142 


0.01565 


0.01311 


0.00931 


0.75 


0.12970 


0.09044 


0.06888 


0.04376 


0.03631 


0.02664 


0.01804 


0.01471 


0.01189 


0.8 


0.11877 


0.09250 


0.06426 


0.04987 


0.04066 


0.03313 


0.02590 


0.02382 


0.01723 


0.85 


0.09331 


0.08519 


0.06628 


0.06425 


0.05193 


0.04726 


0.04363 


0.04031 


0.03580 


0.9 


0.12622 


0.11055 


0.09563 


0.08920 


0.08563 


0.08485 


0.08410 


0.08413 


0.08612 


0.95 


0.12438 


0.11344 


0.10984 


0.11089 


0.11789 


0.12669 


0.14009 


0.15206 


0.17450 



Table 3: - H\ ■ n°--'{hin)-°-'^ 



n points 



H 


100 


250 


1000 


2500 


10* 


2.5 • 10* 


10^ 


2.5 • 10^ 


10« 


0.55 


0.59403 


0.56033 


0.38778 


0.39788 


0.38848 


0.41418 


0.27528 


0.34462 


0.32254 


0.6 


0.51028 


0.42327 


0.40402 


0.34882 


0.34127 


0.35107 


0.30255 


0.29419 


0.28218 


0.65 


0.54879 


0.49122 


0.36011 


0.35142 


0.31289 


0.31248 


0.29686 


0.30960 


0.27064 


0.7 


0.37926 


0.34744 


0.36932 


0.33076 


0.30417 


0.25725 


0.25904 


0.26879 


0.26421 


0.75 


0.45857 


0.38776 


0.39935 


0.31116 


0.35343 


0.31994 


0.29864 


0.30156 


0.33742 


0.8 


0.41990 


0.39658 


0.37253 


0.35462 


0.39573 


0.39783 


0.42871 


0.48825 


0.48905 


0.85 


0.32990 


0.36523 


0.38428 


0.45683 


0.50541 


0.56751 


0.72214 


0.82624 


1.01618 


0.9 


0.44627 


0.47396 


0.55444 


0.63424 


0.83341 


1.01898 


1.39207 


1.72453 


2.44410 


0.95 


0.43976 


0.48637 


0.63677 


0.78848 


1.14744 


1.52138 


2.31880 


3.11678 


4.95257 
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